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The geometry around a rotating massive body, which carries charge and electrical currents could 
be described by its multipole moments (mass moments, mass-current moments, electric moments, 
and magnetic moments). When a small body is orbiting around such a massive body, it will move 
on geodesies, at least for a time interval that is short with respect to the characteristic time of 
the binary due to gravitational radiation. By monitoring the waves emitted by the small body we 
are actually tracing the geometry of the central object, and hence, in principle, we can infer all its 
multipole moments. This paper is a generalization of previous similar results by Ryan. In his paper 
Ryan explored the mass and mass-current moments of a stationary, axially symmetric, and reflection 
symmetric, with respect to its equatorial plane, metric, by analyzing the gravitational waves emitted 
from a test body which is orbiting around the central body in nearly circular equatorial orbits. In our 
study we suppose that the gravitating source is endowed with intense electromagnetic field as well. 
Due to its axisymmetry the source is characterized now by four families of scalar multipole moments: 
its mass moments Mi, its mass-current moments Si, its electrical moments Ei and its magnetic 
moments Hi, where I = 0, 1,2, . . .. Four measurable quantities, the energy emitted by gravitational 
waves per logarithmic interval of frequency, the precession of the periastron, the precession of the 
orbital plane, and the number of cycles emitted per logarithmic interval of frequency, are presented 
as power series of the newtonian orbital velocity of the test body. The power series coefficients are 
simple polynomials of the various moments. If any of these quantities are measured with sufficiently 
high accuracy, the lowest moments, including the electromagnetic ones, could be inferred and thus we 
could get valuable information about the internal structure of the compact massive body. The fact 
that the electromagnetic moments of spacetime can be measured demonstrates that one can obtain 
information about the electromagnetic field purely from gravitational wave analysis. Additionally, 
these measurements could be used as a test of the no-hair theorem for black holes. 
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I. INTRODUCTION 

While, it is well known for quite some time that the 
geometry of the vacuum around a massive object is re- 
lated directly to the multipole moments 0, 0, SE| of the 
central object, as in newtonian gravity, Ryan |5| was the 
first who attempted to map the spacetime geometry of a 
central axisymmetric body, through its mass and mass- 
current moments, on a few measurable physical quanti- 
ties that are related to the kinematics of a hypothetical 
test body that is orbiting around the central object while 
emitting gravitational waves. Later, Ryan |(| used the 
outcome of his work to perform the analysis on the out- 
put of gravitational wave detectors, in order to extract 
the moments of the central body around which a much 
lower-mass body is orbiting and emitting waves. 

Fortunately, we live at an era where technology may 
give us the opportunity to observe astrophysical phenom- 
ena that are related to the highly distorted geometry of 
a compact massive central object. Especially, the detec- 
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tion of gravitational waves from binaries will provide us 
with data that are strongly dependent on the geometry 
itself. Apart from the Earth-based detectors which are 
restricted to detect moderate mass binaries (e.g., 1 to 
300 M Q for LIGO) 7], the Laser Interferometer Space 
Antenna, LISA, when it will fly in space and start oper- 
ating, is expected to explore the geometry of very mas- 
sive objects with exceptional high precision Besides, 
telescopes with higher and higher resolution, operating 
at various regions of the electromagnetic spectrum, are 
able to probe the close neighborhood of compact massive 
objects, as in accretion disks and jets of AGN's. All this 
information which is related, by one way or another, to 
the motion of small objects in the curved geometry of 
massive astrophysical objects could somehow shed light 
on the internal structure of the central object. Of course, 
it is not expected to fully determine its structure from 
the knowledge of its moments (this is not possible even 
in Newtonian gravity). However, knowing the moments 
of the central body, could set restrictions on the various 
models that are assumed to describe the interior of the 
central body. Moreover, the no-hair theorem in the case 
of a large black hole at the role of the central object could 
be fully tested, if from our moment-extraction formulae 
we get the values of electric and magnetic moments. 
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The binaries we have considered in our paper are ideal- 
ized with regard to the following quite realistic assump- 
tions: 

(i) The central object is assumed to be stationary and 
axisymmctric, and characterized by reflection symmetry 
with respect to its equatorial plane. This is expected 
to be true for a quiescent massive object with inter- 
nal fluid motions that are strictly toroidal. The ax- 
isymmetry gives as the freedom to describe the space- 
time geometry with scalar, instead of tensorial, multi- 
pole moments (see Ref. 0). The same is true in New- 
tonian gravity as well, although in that case all other 
moments except the mass moments do not show up 
in the expansion of the gravitational potential. Since 
we are taking into account the electromagnetic content 
of the central object as well, we are considering four 
families of multipole moments to characterize the ge- 
ometry around the central body: Its mass moments 
Mo, M%, M4, . . ., its mass-current moments Si, S3, . . ., its 
electric charge moments Eo, Ei, E%, . . ., and its magnetic 
moments Hq , Hi , H2 , Especially, the Mq = M mo- 
ment is the mass of the object, Si is its angular mo- 
mentum, Eq = E is its charge, and H\ is its magnetic 
dipole moment. In every family of moments, each mo- 
ment appears in steps of 2, and this holds good for the 
electromagnetic moments as well Q. This property is 
due to the reflection symmetry of the metric itself (cf., 
Q). The geometry of such objects can be described by 
the Papapetrou metric which consists of only two dy- 
namical functions. The third one (see Section |HJ could 
be easily inferred from the first two. 

(ii) Although we plan to extend our exploration in a 
generic geodesic motion around such a central object, in 
the present paper we only take into consideration nearly 
equatorial and nearly circular geodesic orbits of test bod- 
ies in the fixed geometry of a central massive object. We 
know that gravitational radiation from a test body that 
is far from its innermost stable circular orbit tends to cir- 
cularize the orbit [ToL ITl| , and therefore the orbit could 
safely be considered circular if it has sufficiently long time 
to evolve without being perturbed by other objects. Also, 
we know that at least for not extremely fast rotating Kerr 
black holes the evolution of non-equatorial orbits due to 
radiation reaction is such that their inclination remains 
almost constant while their radius decreases [T^|. 

(iii) The energy emitted in the form of gravitational 
radiation will be assumed to be given by the quadrupole 
formula, since there is no known way to fully analyze 
the wave emission in a generic geometrical background. 
Furthermore, we assume that this energy is carried away 
by waves at infinity, and there is no energy loss through 
any horizon, or due to thermal heating of the surface of 
the central object from the impact of gravitational waves. 

The rest of the paper is organized as follows: First, 
in Sec. |n] we define the observable quantities that will 
be used to measure the moments of the metric. These 
quantities are the periastron precession fi p , the preces- 
sion of the orbital plane Q z , the energy emitted at in- 



finity per logarithmic frequency change AE/fi, and the 
number of cycles of the primary gravitational waves per 
logarithmic frequency change AA. Especially the latter 
one, which can be measured with high accuracy by the 
broad-band wave detectors that are operating now, or 
will be built in the near future, is computed assuming 
that the phase of the waves is coming simply from the 
dominant frequency, / = Q/n. We also show how one 
defines the mass and electromagnetic multipole moments 
of an axially symmetric body with reflectional symmetry, 
and how these moments uniquely determine the metric 
of the space around the object. For the equatorial plane, 
and slightly out of this plane, we write the metric as a 
power series of the Weyl radial coordinate p. The coeffi- 
cients of the power series are polynomials of the various 
moments of spacetime. We end up this section by ex- 
plaining the implications of reflection symmetry of the 
metric on the electromagnetic fields. Having in hand all 
these expressions that connect the observable quantities 
to the metric and hence to the moments, we proceed 
in Sec. II I II to write down expressions for the four as- 
trophysical quantities, as power series of v = (Mfi) 1 / 3 , 
where M is the mass of the central object and is the 
orbital frequency of the test body, observed at infinity. 
This quantity is the newtonian orbital velocity of the or- 
biting body, and is a measure of the gravitational field 
strength. Following the analysis of Ryan [5| , we present 
the power expansion of AA with coefficients that include 
only the leading order contribution of each central-body 
multipole moment. It is argued that these are the only 
terms we could get without getting into the complicate 
analysis of wave emission. Finally, in Sec. II VI we comment 
on how the gravitational wave analysis could inform us 
about the moments of the central object, the accuracy 
with which these moments could be computed, and the 
implications they could have on the observational verifi- 
cation (or not) of the full no-hair theorem (when charges 
are included) . Throughout the paper units are chosen so 
that G = c = 1. 



II. OBSERVABLE QUANTITIES AND 
MOMENTS 



In this section we briefly present and discuss the mea- 
surable quantities that Ryan has used as tracers to mea- 
sure moments. A thorough presentation and analysis of 
them can be found in Ref. [5|- Then, we present all formu- 
lae that determine the various moments and relate them 
with the metric describing the geometry around a central 
object. Finally, we discuss what kind of electromagnetic 
fields are consistent with the symmetries assumed for the 
metric. 
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A. Quantities that can be measured through 
gravitational-wave analysis 

As is explained in |5j there are four physical quantities 
in a binary with high-mass-ratio, that can, in principle, 
be measured through the gravitational radiation emitted 
by the binary, and are straightforwardly related to the 
spacetime metric of the massive body. It is exactly these 
quantities that we will use here as a basic information 
in order to extract the various multipole moments of the 
central object. These are: (i) Q p , the periastron pre- 
cession of the low- mass body, (ii) Vt z , the orbital-plane 
precession of the low-mass body, (hi) AE, the energy 
emitted as gravitational waves per logarithmic interval 
of frequency, and most important (iv) AN, the num- 
ber of gravitational wave cycles per logarithmic interval 
of frequency. The first two quantities are computed by 
analyzing the geodesic motion of nearly circular, nearly 
equatorial orbits of a test body on a fixed background 
metric, and can be measured through the modulation 
they induce on the gravitational waves emitted by the 
binary. The third one, although it is directly related to 
gravitational radiation emitted by the binary, it can be 
easily inferred by the functional relation of the energy 
of the test body, which changes adiabatically, to its or- 
bital frequency, which is simply half the primary wave 
frequency. Finally, the fourth one is the best measurable 
quantity, since the phase matching used in data analy- 
sis leads to a highly accurate estimation of the frequency 
dependence of phase, assuming detection has been es- 
tablished. The computation of this quantity involves a 
number of approximations, since it is directly related to 
the mechanism of gravitational wave emission on a com- 
plicated metric background. 

Here, for the sake of completeness, we rewrite the ex- 
pressions of Ryan [Bj], that relate all these observable 
quantities with the metric without further comments on 
how these are computed. Both precession frequencies are 
given by 
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where a stands for p, or z. Actually, the frequencies 
written above correspond to the difference between the 
orbital frequency and the frequency of perturbations in 
p, or z, since these differences are expected to show up in 
gravitational waves as a modulating frequency. Of course 
these frequencies are accurate only for orbits that are 
slightly non circular and slightly non equatorial, other- 
wise the frequencies would depend not only on the metric 
but on specific characteristics of the orbit, its eccentricity 
and inclination. 



The energy per unit test-body mass for an equatorial 
circular orbit in an axially symmetric spacetime is 
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and thus, the specific energy released as gravitational 
radiation per logarithmic interval of frequency is 
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The expression above assumes that all the energy lost 
from the test body has been emitted at infinity as ra- 
diation, and is neither wasted as thermal energy on the 
fluid of the central object, nor is it "lost" through any 
horizon. 

The number of gravitational-wave cycles spent in a log- 
arithmic interval of frequency is 
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where dE wstvc /dt is the gravitational- wave luminosity, 
which is assumed to be exactly the rate of energy loss of 
the orbiting test body. As Ryan has analytically shown 
the main contribution of dE wavc /dt comes from the mass 
quadrupole radiative moment of the binary. More specif- 
ically, up to fourth order of v = (Affi) 1 / 3 after the lead- 
ing order, the gravitational wave luminosity is accurately 
computed from the quadrupole formula 
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plus a contribution of the current quadrupole radiative 
moment, due to the motion of the central object around 
the center of mass, 
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see Ref. We should also put by hand further addi- 
tional contributions of dE w&vc /dt, due to post-newtonian 
corrections. The corresponding contributions, up to v A , 
are simply numerical if one computes them from pertur- 
bation analysis in a Schwarzschild background. A com- 
parison though, between the final formula (55) of for 
dE/dt, and formula (3.13) of 13], which is based on per- 
turbations on a Kerr background, shows that at least up 
to v A order, the terms of |5| that include S\ and M%, 
which come from the corresponding contributions that 
are given by Eqs. (|5I6() . in the case of a Kerr metric, 
are equal to the ones of Ref. ^| . This agreement indi- 
cates that up to v A order we can simply add the follow- 
ing numerical post-newtonian terms (the corresponding 
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terms of if we set q = 0) in the rest contributions of 
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Finally, in order to compute the number of cycles AN, 
one has to add up all contributions of dE wavc /dt and 
combine them with the expression for AE/pi, which is 
given above (see Eq. (J3J). 

To get expressions for all these quantities, that are 
straightforwardly connected to the moments of the cen- 
tral object, and can be observationally measured, first 
one has to reexpress the metric functions in terms of all 
moments. Also the relation between the radius p and the 
orbital frequency of the test body SI, through moments, 
is necessary so as to finally express the four measurable 
quantities as power series of v. 



where V denotes the gradient in a cartesian 3D space 
(p, z, 4>) and 5R(. . .), S(. . .), here and henceforth, denote 
the real and imaginary part, respectively, of the complex 
function in parentheses. An asterisc * denotes complex 
conjugate. The third metric function, 7, is then easily 
computed by integrating the partial derivatives d"//dp, 
d"//dz, which are given as functions of derivatives of all 
other functions 
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B. Moments describing spacetime 

In our paper, we consider only stationary axisymmetric 
objects that are symmetric with respect to their equato- 
rial plane. These symmetries are more or less realistic 
assumptions for a quiescent massive rotating astrophys- 
ical body around which much smaller bodies orbit. The 
metric of such a central object alone could be written in 
{Up, z, 4>) coordinates, in the form of Papapetrou metric 



dz 2 



ds 2 = -F(dt - uj d<p) 2 + 

1 1 

where F,lu, and 7 are the three functions that fully de- 
termine a specific metric. These are functions of p and 
\z\ only, due to axisymmetry and reflection symmetry. 
Einstein's equations in vacuum guarantee that once F 
and uj are given, 7 can be easily computed (see (15^. 
Once we incorporate an electromagnetic field in the vac- 
uum around the compact object, the source of which 
is the compact object itself, which allows spacetime to 
have the same symmetries, the metric above still de- 
scribes the electrovacuum spacetime, but now the metric 
and the electromagnetic field should satisfy the Einstein- 
Maxwell equations. In order to fully compute the metric 
functions one more complex function, <f>, which is related 
to the electromagnetic field, is necessary. F, uj, and <& 
themselves can be determined by solving the so-called 
Ernst equations [l6l Il7| . which are nothing more than 
the Einstein-Maxwell equations in a different form. It is 
a system of non-linear complex differential equations of 
second order: 
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The relation of the two complex functions, £ and $, 
with the metric functions is the following 
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where tp is related with g t( p through 
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Note that there is a sign difference in Eq. (22) of 
which has been corrected in a later paper of Ryan [li 
and comes from an odd convention of uj used by Ernst 
(see relevant comment of [T^). 

Instead of £ and one could use two new complex 
functions £ and q, that play the role of gravitational po- 
tential and Coulomb potential respectively, and are more 
directly connected to the mass and electromagnetic mo- 
ments of the central body. These potentials are related 
to the Ernst functions by 
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and can be written as power series expansions at infinity 
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where 



where are computed recursively by 



p 2 + 



(18) 



and a.ij , bij are coefficients that vanish when i is odd. 
This reflects the analyticity of the potentials on the z- 
axis. The tilded quantities, here and henceforth, are the 
conformally transformed ones, which are essential for cal- 
culating the moments (see []}). 

Due to Ernst equations I|9I10|I the above power expan- 
sion coefficients, <Zjj and bij, are interrelated through the 
following complicated recursive relations 
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where m = r — k— p,0<k<r,0<p<r — k, with 
/c and p even, and n = s — I — g, < Z < s + 1 , and 
— 1 < g < s — I. Essentially, these relations are simply 
an algebraic version of Einstein-Maxwell equations for 
the coefficients of the power expansion of the metric and 
the electromagnetic field tensor. The recursive relations 
(|19I20|1 could be used to build the whole power series of 
£ and q from a full knowledge of the metric on the axis 
of symmetry 
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In [2(j a method of calculating the complex multipole 
moments of the central object in terms of the m^'s and 
qiS is presented. In brief, the gravitational moments are 
given by 
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The electromagnetic moments Q n are computed from 
exactly the same formulae, by simply replacing the initial 
term with q instead of £. flu, R12 and i?22 are given 
by 
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and 71 = y p and 72 = j, z can be expressed in term of 
Rij as 
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~p(R n - R22) 



(26) 



The mass moments M n and the mass-current moments 
S n are related to P n by 



P n — M n + iS n 



(27) 



whereas the electric moments E n and the magnetic mo- 
ments H n are related to Q n by 

(28) 
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Since this algorithm can be used to evaluate the mo- 
ments in terms of the m^'s and q^s, one can invert these 
relations and express the m,'s and q^s in terms of the 
moments: 

m n = a 0n = M n + iS n + LOM, 
q n = b 0n = E n + iH n + LOM, (29) 

where "LOM" stands for lower order multipole moments 
of any type. Thus, we can use the recursive relations i|19|) 
and l|20[) to evaluate the a,j and bij coefficients in terms of 
the moments. Finally, following the procedure presented 
in the beginning of this subsection (Eqs. I|11I17|I % ) we can 
express the metric functions, and their first and second 
derivatives as power series of p and z with coefficients 
that are simple algebraic functions of the moments of 
the massive body. Since in our study we have confined 
the motion of the test particle on the equatorial plane, 
we actually need to compute everything at z — which 
makes calculations far simpler than what they seem. 

C. Reflection symmetry and electromagnetic 
moments 

Following Ernst 0, and using Papapetrou's met- 
ric © we end up with the following Einstein-Maxwell 
equations: 

V ■ [p- 2 F (VA 3 - ojVA 4 )] = 0, (30) 
V • [F _1 VA 4 + p - 2 Flu (VA 3 - uVA A )} = 0, (31) 
V • [ p - 2 F 2 Vuj - 4p- 2 FA 4 (VA 3 - ujVA 4 )] = 0, (32) 
FV 2 F = VF • VF - p-' 2 f 4 Vuj ■ Vlu + 2FVA t ■ VA 4 + 
+ 2p 2 F 3 (VA 3 - uVA A ) ■ (VA 3 - uVA A ) , (33) 

where A 3 and A A denote the A$ and A t components of 
the electromagnetic 4-potential respectively, and V is the 
three dimensional divergence operator in Weyl coordi- 
nates. As we have already mentioned, the three metric 
functions F, uj and 7 are functions of p and \z\ only, due 
to the assumed symmetries. From the equations above 
one cannot easily tell what are the symmetries, if any, 
that are inherited in A3 and A A . It is obvious though 
that A3 and A4 being either both odd or both even func- 
tions of z is consistent with the reflection symmetry of 
the metric functions. We will argue that these are the 
only reasonable types of the electromagnetic field. 

The Ernst potential $ is defined by Ernst through the 
following relations: 

p- 1 f(VA 3 -iuVA 4 )^nxVA / 3 , (34) 



$ = A 4 + iA 3 . (35) 

where n is the unit vector in the azimouthal direction. 
From Eq. I|14|). we obtain the metric function g tr j,. Since 



we want to end up with metric functions that are even 
functions of z, then the whole intergrand should be even 
as well. This could be acomplished if p is an odd func- 
tion of z and the real and imaginary part of $ are cither 
even and odd or odd and even functions of z respectively. 
If none of the above holds then in order to get an even 
function for g t( j >1 one must impose a restraining functional 
relation between tp and But tp and should be in- 
dependent in order to describe a generic spacetime with 
the symmetries mentioned above. Therefore, by virtue 
of Eqs. I|34I35|) . either A4 and A 3 are both even, or both 
odd functions of z. Consequently, for both cases the 9?(£) 
is an even and the 3(£) is an odd function of z. 

Now, we can use this information to conclude that, 
by virtue of Eq. I|15|) . the action of reflection symme- 
try leaves invariant and changes the sign of 
Thus only even order mass moment and odd order cur- 
rent mass moments will occur [2lJ. For the electromag- 
netic moments things are not univocal since we have two 
discernible electromagnetic cases that are consistent with 
the symmetries of the metric. If both A 3 and A4 are odd 
functions of z then the action of the reflection symme- 
try leaves 3?(g) invariant and reverses the sign of Q(q), 
which means that only even order electric field moments 
and odd order magnetic field moments will occur. On the 
other hand if both A3 and A A are even functions of z then 
the action of the reflection symmetry leaves $s(q) invari- 
ant, but reverses the sign of 3t(q), which means that in 
this case only odd order electric field moments and even 
order magnetic field moments will occur. 



III. THE POWER EXPANSION FORMULAE 

Combining all formulae that are given in Sec. [H] even- 
tually we can express all four measurable quantities as 
power series of v = (Mfi) 1 / 3 with coefficients that have 
explicit dependence on all four types of moments. The 
choice of v as a dimensionless parameter to expand all 
physical quantities is warranted from the fact that the 
inspiral phase of a binary, the best exploitable part in 
gravitational- wave analysis p^. involves comparatively 
low magnitudes of v. All measurable quantities have been 
transformed to a dimensionless form as well, for example 
by dividing the two frequencies Q p , Q z , by the orbital 
frequency £1. 

Since the metric functions and their derivatives, are 
expressed as functions only of p at the equatorial plane, 
in order to express all measurable quantities as power 
series of v, we need also a power series expansion of p 
with respect to v, or equivalently Q. Thus we have to 
invert the function f2(p), at least as a power expansion. 
From an elementary analysis of circular geodesies on the 
equatorial plane (see 0) we know that 



-ff«0, P + y/(.9t0, P ) 2 - {gtt,p)(g^, P ) 
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In the following part of this section we explain the al- 
gorithm that one should follow, in order to obtain the 
power series for f2 p /f2, Vt z /Vt : AE/p, and AN. One 
starts with a power series of £ and q of the form given 
by Eq. I|17|l . Since no higher than second derivatives 
of the metric functions with respect to z are necessary, 
one should keep ay's and fry's with < j < 2, and as 
many values of i as one needs to carry the power series 
expansion of the measurable quantities at a desirable or- 
der. In our paper where all quantities are written up 
to no higher than v 11 order, we only need i's in the in- 
terval < i < 4. All quantities that are expressed as 
power series of p and z, are evaluated at z — z = at 
the end, and thus, all expressions are finally power series 
of p = 1/p, due to Eq. (|18[1. Although the ay and fry- 
are polynomials of various moments, from the practical 
point of view it is preferable to keep them as they are, 
and replace them by their moments dependence only at 
the final expressions. Then from £ and q we construct £ , 
and F, ip (cf., Eqs. jl3H5jl ). These are sufficient to 
build all metric functions through Eqs. I|11II2II4I1 . Next, 
following the procedure described above, we expand Q as 
a power series of 1/p, by virtue of Eq. (JJSJ). This series is 
inverted and in this way we obtain 1/p as a power series 
of O, which then can easily be turned into a power series 
of the dimensionless parameter v. 

Now, the power series representing 1/p will replace all 
1/p terms appearing at the expansions of the metric, its 
derivatives, and all other physical quantities depending 
on them (Eqs. Illl5fl '). Finally, one has to rewrite the ay 
and fry terms appearing at the coefficients of all these 
power series as polynomials of the various moments. The 
recursive relations (|I9I2U[1 relate all ay and fry with = 
aofc and qk = fro/c , which are directly related to the scalar 
moments of spacetime through Eqs. (24,25) of Ref. pcj . 

The algorithm described in the previous two para- 
graphs has been carried out with Mathematica, and has 
been checked for the following two subcases: (i) When all 
electromagnetic fields are turned off, by erasing all elec- 
tromagnetic moments (Ei = Hi = 0), our expressions for 
O p , Vtz,AE/p, AN are identical to the ones computed by 
Ryan [f| . (ii) For the Kerr-Newman metric it is quite easy 
to compute O, tt p , Cl z , and AE / p for a quasi-equatorial, 
quasi-circular orbit. Actually, there is no need to use 
Weyl coordinates to describe the metric; one could sim- 
ply work with the metric in the usual Boyer-Lindquist 
coordinates (see Eq. (33.2) of 23]), and compute ev- 
erything according to the formulae given above, by re- 
placing the derivatives with respect to z, with the corre- 
sponding derivatives with respect to 9 around 9 = tt/2. 
The expressions for p though should be replaced with 
iSts ~ Qttg^xp) 1 ! 2 ■ Thus, if in the power series for Q p , O z , 



and AE/p that are written below (cf., Eqs. 1 (351 139140(1 ). 

one makes the following substitutions for the moments 

M 2l = (-l) l Ma 21 , S 2l+1 = (-l) L Ma 2l +\ 

E 2 i = (Q/M)M 2h H 2l+1 = {Q/M)S 2l+1 , (37) 

according to pp|. the expressions we obtain are identi- 
cal to the ones obtained directly from the Kerr-Newman 
metric. 

There is one more thing that should be pointed out 
before we write down the power series for all four ob- 
servable quantities. As is explained in Sec. Ill ( 'I there are 
two possible cases for the electromagnetic field that lead 
to reflection-symmetric spacetimes. The first case (with 
odd A3 and A4 as functions of z) is the one that de- 
scribes an electric field that is reflection symmetric (like 
in a monopole electric field), and a magnetic field that is 
reflection antisymmetric (like in a magnetic dipole field). 
Henceforth we shall call this case electric-symmetric case 
(es). In that case only the even electric moments and 
the odd magnetic moments show up in the moment anal- 
ysis of spacetime. Thus fro; is real for even I's and purely 
imaginary for odd i's. The other case (with even A3 and 
A4 as functions of z) is the one that describes an elec- 
tric field that is reflection antisymmetric (like in a dipole 
electric field), and a magnetic field that is reflection sym- 
metric (like in a magnetic quadrupole field). Henceforth 
we shall call this case magnetic-symmetric case (ms). In 
that case only the odd electric moments and the even 
magnetic moments show up in the moment analysis of 
spacetime. Thus fro; is real for odd I's and purely imagi- 
nary for even i's. Although classically we do not expect 
the central object to carry any magnetic monopole, the 
zcroth order magnetic moment shows up formally in the 
terms of a generic magnetic-symmetric case, and thus we 
have not omitted it. 

The power series expansion for Q p , Cl z , and AE/ p take 
the following form: 

Q °° 

= £fl»« n . (38) 

n=2 

o 00 

^ = £^"' (39) 

n— 3 

— = f>< ( 4 °) 

r n=2 

while the corresponding coefficients, up to 9nth order for 
the two frequencies and up to Ilnth order for the radiated 
energy, in the two distinct electromagnetic cases (es) and 
(ms) are 
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We note that the higher order terms we have computed 
are by one order lower than the corresponding higher 
order terms of Ryan. This is due to computing power 
restrictions, since as one goes to higher order terms our 
coefficients become far richer in moments than Ryan's. 
The fact that we have two new sets of moments (the 
electromagnetic ones) allows many more combinations of 
moments in high order terms. Actually, from a practical 
point of view these expansions are far more advanced 
than what will be used in gravitational wave data analysis 
in the near future. On the other hand the expressions 
above present an important feature: in every new order 
term a new moment shows up. This suggests that a very 
accurate observational estimation of the series could in 
principle reveal any moment. 

A glance at the corresponding terms of the two electro- 
magnetic cases shows that each combination of moments 
for the (es) case is numerically equal to the corresponding 
combination for the (ms) case, if the electric and mag- 
netic moments are interchanged. The sign though is the 
same for combinations of pure electric, or pure magnetic 
moments, but opposite for combinations of electric and 
magnetic moments. 

The difference by two in the order of the highest com- 
puted order term between the power series for the O's 
and AE is due to the second derivatives that appear in 
Eq. p. 

Finally, in order to express AN also as a power series of 
v, we need to expand dE wavc /dt as power series of v. As 
was explained in Sec. Ill Al we cannot work out the pertur- 
bative analysis of gravitational wave emission at a generic 
spacetime background; we can only obtain accurate ex- 
pressions for dE wave /dt up to v 4 after the leading order. 
However, we know the numerical factor of the higher 
order moment appearing at any higher than v 4 order. 
These higher order moments come from the power series 
expansion of p itself through Eq. JSJ which describes the 
main contribution to energy radiation. All other contri- 
butions depend on lower moments at the same order of v. 
Therefore, in the following formulae for p and dE wavc /dt, 
we write the power series coefficients explicitly up to the 
fourth order, while instead of giving the explicit form 
of all higher order coefficients, we give only the higher 
moment term that occurs at each order of the power ex- 
pansion. More specifically, in order to compute the power 
expansion of dE wajVC /dt we add up all three power series 
contributions of Eqs. (|5I6I7[1 . Thus, we yield 



P = Mv- 2 (l + J2p n v n 



(41) 



and 



dE wavc 32 / fi 



(&) » 10 KEv 



(42) 



where the p n and W n coefficients for the two electromag- 



netic cases are respectively 
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(43) 
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(46) 



In the expressions above the indices k run from 1 to infin- 
ity, and the term LOM is an abbreviation for lower order 
moments that appear at a specific order in the expan- 
sion. We note that in all these coefficients the same fea- 
ture with respect to the corresponding terms in the two 
electromagnetic cases that was mentioned before arises. 

By combining the power series for AE//I (Eq. ^Ai)\i ) 
with the one for dE wavc /dt (Eq. (|42|l ) we obtain the power 
series expansion of AN: 
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(47) 



where the N n coefficients for the two electromagnetic 
cases are given by the following polynomials of the mo- 
ments 
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As in the rest three measurable quantities, the power 
expansion of A./V is such that in every order term a new 
moment, which was not present in any lower order term, 
occurs. This proves that all moments can in principle be 
unambiguously extracted from accurate measurements of 
AN. 



IV. USING THE RESULTS IN 
GRAVITATIONAL WAVE ANALYSIS 

Although wc have not quantitatively explored the im- 
plications of our results on the estimation of errors in 
determining the various moments from a gravitational- 
wave data analysis, as it has been done by Ryan in 0, 
we could make some general comments. Actually, the 
only difference of our results from the ones of |5| is that 
more moments are showing up at each coefficient in the 
power expansions of all observable quantities, and thus 
Ryan's estimates for each term apply equally well here. 

As is shown in 6] the first generation of LIGO is not 
expected to be able to extract the first two moments (Si 
and M 2 ) with high accuracy (~ 0.05 for the former and 
~ 0.5 for the latter one), by analyzing the phase of the 
waves. If we allow for electromagnetic fields as well, the 
corresponding monopole (which classically is expected to 
be very close to zero) will be measured with even higher 
accuracy than the other two mass moments, since the 
charge of the source (or the magnetic monopole in case of 
some exotic body) is present at even lower order, namely 
in the v 2 term, while the electric dipole, or the mag- 
netic dipole, that first show up at the v 5 term will be 
measured with rather disappointing accuracy. On the 
other hand, analyzing the data of LISA leads to accu- 
racies almost two orders of magnitude higher than the 
corresponding for LIGO. Thus, it seems quite promising 
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that LISA will give us the opportunity to measure the 
first few moments, including the electromagnetic dipole 
moments, quite accurately. Also, the fact that in every 
new order in the power series of AN a new moment ap- 
pears is significant, since this means that in principle a 
unique set of moments arises from an accurate estimation 
of all power series terms. Actually, there are two possi- 
ble sets of moments; one for each electromagnetic case, 
since we cannot a priori exclude one of them. We can 
only exclude one of the two sets on physical grounds, if 
only one of them leads to a physically reasonable classical 
object (for example a highly magnetized compact object 
is physically preferable to a compact object with a huge 
electric dipole). If we manage to measure a few lower 
moments, we can check if they are interrelated as in a 
Kerr-Newman metric po| . A positive outcome of such 
a test will be of support to the black-hole no-hair theo- 
rem in the case that the central object is a black-hole. 
The case of observational violation of the black-hole no- 
hair theorem could cither mean that the central object is 
not a black-hole, or that the theorem does not hold. Of 
course to assume the latter an extra verification that the 
central object is indeed a black hole is necessary. Indi- 
cations that the central massive compact object is not a 
black hole would imply the existence of an exotic object 
(e.g., soliton star, naked singularity, etc.). If the central 
body's mass is measured to be within the stellar limits 
(e.g., a massive neutron star) we could get valuable infor- 
mation about its electromagnetic field, like its magnetic 
dipole field. 

While the phase of a gravitational wave is the quantity 
that can be most accurately measured, since a large num- 
ber of cycles (a few thousand for LIGO and a few hundred 
thousand for LISA in case of binaries with high-ratio of 
masses) is sweeping up the sensitive part of the detectors, 
the two precession frequencies tt p and Q z , can in princi- 
ple be measured if the detectors become more sensitive 



and templates that describe modulating waves are used 
|24|. If this ever become possible one could use any of 
them to test the no-hair theorem. This would demand 
no more than the four lower order terms, since accord- 
ing to this theorem all moments depend on only three 
quantities (mass, angular momentum, and total charge). 
Actually, from measurements of modulating frequencies 
we could not at first determine which frequency corre- 
sponds to each precession. However, the power expan- 
sions of the two frequencies begin at a different order, 
and thus we could discern them. Unfortunately, every 
term of these expansions contains more than one first oc- 
curring moment. However, expansions of f2 p and Q z , if 
used simultaneously, could lead to the full determination 
of the moments. 

Our analysis demonstrates that it will be possible to 
determine all types of multipole moments of the central 
object, from future gravitational wave measurements. 
Thus, apart from spacetime geometry, we could also de- 
termine the central body's electromagnetic fields. Al- 
though the data of LISA should be suitable for extracting 
such information with high accuracy, the assumptions of 
circular and equatorial orbit are not that realistic. From 
this point of view we consider our work as a step towards 
a more detailed analysis with not so restrictive assump- 
tions. 
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